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Abstract 

We calculate the supersymmetric 0{a s ) QCD corrections to the width of 
the decay H + — > tb within the Minimal Supersymmetric Standard Model. We 
find that the QCD corrections are significant, but that they do not invalidate 
our previous conclusion at tree-level on the dominance of the tb mode in a wide 
parameter region. 



1 Introduction 

In the Minimal Supersymmetric Standard Model (MSSM) [p], @] two Higgs doublets 
are necessary, leading to five physical Higgs bosons h°,H°,A°, and H ± || |[. If 
all supersymmetric (SUSY) particles are very heavy, the charged Higgs boson H + 
decays dominantly into tb; the decays H + — > t + u and/or H + — > W + h° are dominant 
below the tb threshold |3], [§] . In ref . M all decay modes of H + including the SUSY- 
particle modes were studied in detail; it was shown that the SUSY decay modes 
H + — > tibj(i,j = 1,2) can be dominant in a large region of the MSSM parameter 
space due to large t and b quark Yukawa couplings and large t- and 6-mixings, and 
that this could have a decisive impact on H + searches at future colliders. Here ti (bj) 
are the scalar top (scalar bottom) mass eigenstates which are mixtures of ti and t R 
(b L and b R ). 

The standard QCD corrections are very large for the width of H + — > cs and 
can be large (+10% to —50%) for that of H + — > tb 0]. The QCD corrections from 
the SUSY-particle loops are calculated within the MSSM for H + — > tb in and 
turn out to be non-negligible (~ 10%) for certain values of the MSSM parameters. 
This suggests that the QCD corrections to H + — > tb could also be large. Therefore it 
should be examined whether the result in 0] remains valid after including the QCD 
corrections. 

In this paper we calculate the 0(a s ) QCD corrections to the width of H + — > 
t{bj within the MSSM. To the best of our knowledge they are not known in the 
literature. We obtain the complete 0(a s ) corrected width in the DR renormalization 
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scheme (i.e. the MS scheme with dimensional reduction [§) including all quark mass 
terms and q~i — q~R mixings. The main complication here is that the q~L — q~R mixing 
angles are renormalized by the SUSY QCD corrections. We find that the corrections 
to the tb width are significant but that the tb mode is still dominant in a wide 
parameter range. 

2 Tree level result 

We first review the tree level results 0. The squark mass matrix in the basis (qi, 
(Jr), with q = i or b, is given by 0, f| 

2 \ / ™2 



where 



m 2 LL = MQ + m 2 q + m| cos 2(3 (I q - Q q sin 2 6 W ), (2) 
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m RR 



= Mrftfo + m q + m| cos 2(3Q q sin 2 W , (3) 



2 _ 2 / m(A - yucot/3) (q = t) 

U>L!; ' " >liL \ /,„,(. 1,,-,, tan 1) (q b) ' 



(4) 



and 



q ( cos 6q sin 9q \ , , 

ia \ -sm9g cos9 q J ■ { ! 

Here the mass eigenstates q~i(i = 1,2) (with m qi < m q2 ) are related to the SU(2)^ 
eigenstates q a (a = L,R) as % = R q ioi q a - Note that in the sign convention used here 
the parameters A tjb correspond to (— A t ^) of ref. ||. 

The tree-level decay width of H + — ► tj&j is then given by (see Fig. la) 

Y^(H^Ub ] ) = -^ Y \G l] \\ (6) 



where ttih is the H + mass, k = ^mfj, m~ , m| ), y, z) = ((x — y — z) 2 — Ayz) 1 / 2 , 
N c = 3, and 

£ g / m t tan/3 + cot/5 - sin2/3 ra;,(yl 6 tan/3 + //) \ , g.^ 

ij " v^^vk V m t (A iC ot/3 + /i) 2m t m fe /sin2/3 J 1 J 

(7) 

are the H + Ubj couplings || f4j, with g being the SU(2) coupling. 

3 QCD virtual corrections 

The 0(a s ) QCD virtual corrections to H + — > ij&j stem from the diagrams of Fig. lb 
(vertex corrections) and lc (wave function corrections). For simplicity we use in this 
paper the DR renormalization scheme* for all parameters which receive the QCD 
corrections, i.e. m^, A^, and Mq q ^. The renormalized squark mixing angle 9q 



is then defined by the relations (1-5) in terms of the DR parameters m ti &, A t &, and 



The one-loop corrected decay amplitudes G£° rr are expressed as 



G -° rr = G y + «JGg } + tfGg , (8) 

where Gij are defined by (7) in terms of the DR parameters, and SG^ and SG[J^ are 
the vertex and squark wave function corrections, respectively. 

The vertex corrections 8G$ are calculated from the graphs of Fig. lb as 



6G v ~ "17" 



{B {m\, 0, m|) + B (m~ b ., 0, m? ) - B (m#, m| , mf 



* Strictly speaking, our renormalization scheme is the DR scheme |l0(] where the "e-scalar mass" 
is absorbed into M? 



•,U,D 



-2{m H - m~ t . - m- b )C {m- b , A , m f .)}G^ 



-B (m 2 H , m\ , )G fcJ 5j fc 5'J 



+2{(a LL )jj(^t2/2 + mtyi) + {oi RR )ij{mtyi + m b ?/2)}£o(m^, mf, 



+2{(a LL ) i:j m t y 2 + {otLR)ijm g yi + {a RL )ijmgy 2 + (a RR ) i jm t y 1 }Bo(m'l,m 2 g ,m 2 '' 
+2{(a LL ) ij m b y 1 + {otLR)ijm~ g yi + {oiRL)ijm g y 2 + {oiRR)ijm b y 2 }B (ml.,m 2 s ,ml) 
+2{(m t 2 + m 2 - m 2 H )mg((a LR )ijyi + (aju,)^) 
+(m 2 + m~ - ml,)m t ((a LL )ijy 2 + {a RR )ijyi) 
+ (m 2 t +m 2 g- rn\)m b {{a,L L )ijyi + {a RR ) i:j y 2 ) 



+2m g m t m b ({a LR ) ij y 2 + (a^)^^}^^, m|, 



(9) 



where = 4/3, 



/ cos % cos 6*^ 
n 1 1 I — sin % cos 

/ — sin 0£ cos 
aRL ~ ^ -cos0 t -cos0 6 



cos 20g — sin 26*g 
— sin 2^ — cos 2^ 

sin — cos 6i cos ^ 

Gin $~ Gin $~ nnc H~ 



(10) 



cos cos &i — cos 0j sin 9~ b \ f — cos % sin 9 b — ws uus 

sin ^ sin 9 b J ' aLi? I sin 9i sin ^ sin % cos ^ / 

3in%sin0g sin % cos ^ \ , . 
:os^fsin^ cos ^ cos 9~ b J 



— sin 9i cos sin 9\ sin ^ \ / sin 9i e 

^ cos 6>£ cos % sin 0g i ' ~~ y cos 9~ t i 

M — rrib tan (3 = hb sin (3, y 2 = —^2- — m t cot (3 = h t cos /3, (12) 



yi = 

y2m w 



= hb sin (3, y 2 = —jM- — m t cot (3 
\j2m w 



A gluon mass A is introduced to regi 



and irtg is the gluino mass. 

divergences. Here we define the functions A, B , B 1 and C as in 
D) — "fE + log47r): 



;ularize the infrared 
Q (A = 2/(4 - 



A(m 2 



B (k 2 \m\,m 2 ^ 



d D q 



in 2 q 2 
d D q 



77T 



m 2 (A + log(Q 2 /m 2 ) + 1), 



1 



(13) 



m 2 (q 2 — m 2 )((q + k) 2 — m 2 ) 



5 



B x {k 2 \m\,ml) 



C (ml,ml,ml) 



A- f 1 dz]og ^-^l + zr4-z(l-z)^-i6 
Jo Q 2 

ID, 



1 rd u q 



J m 2 (q 2 — m 2 )((q + k) 2 — m 2 .) 
A r 1 , , (1 — z)m 2 + zml — z(l — z)k 2 — i5 . . 

"2 + L dZZl0g ~ — " ' (15) 

d D q 1 
in 2 (q 2 — m 2 )((q + /C2) 2 — Tn?>)((q + p) 2 — m 2 ) 

dx I dy I dz5(l — x — y — z) x 
10 Jo Jo 

{xm\ + ym\ + zm\ — xymf, — yzm 2 — xzm 2 H — id)^ 1 . (16) 



Here p and k 2 are respectively the external momenta of H + and bj, and Q is the DR 



renormalization scale. Note that A is omitted in the DR scheme. 

The squark wave function corrections 5G^ are expressed as 



5G 



i&(m? ) 



11 9 n u t'1 9 

U til bj b j 



2 ^ij' j 



(17) 



where i ^ i' and j 7^ j'. IT?- (A; 2 ) are the one-loop corrections to the two-point functions 
of q^j, which are obtained from the graphs of Fig. lc. Il(fc 2 ) denotes the derivative 
with respect to k 2 . The last two terms in (|T7|) represent the corrections due to the 
renormalization of the g-mixings. The explicit forms are 

a s C F 



47T 



-3£ (m| , 0, m|) - 25 1 (m|, 0, m|) - 4m|S (m|, A 2 



Tra- 



gi' 



-2m|5 1 (m|, 0, m 2 -.) - 4to?5 (to|, m|, m 2 ) - AB^ml , mj, m 2 ) 



-4m?5i(m? ,m?,m 2 ) + 



u-i 



4 sin 29qm q mgBo(m~. , mj, m ? , 



(18) 



and 



47T 



sin46>g(A(m? 2 ) - A(m|J) + 4cos29qm q m g B (m 2 q.,m 2 g ,m 2 q ) 



(19) 



The one-loop corrected decay width in the DR scheme is then given by 

T(H + -> fa) = [|G y | 2 + 2G tJ Re(6G^ + 6G$>)]. (20) 

107T77T< 

Here K po i e refers to k in (6) evaluated with pole squark masses. The width of (20) is 
infrared divergent. 

In the numerical analysis we take the pole quark masses as inputs. The DR 
quark masses are obtained from the pole quark masses by using 

m q(Q)uR = rn g (po\e) ^— ^ [2m q (B (m 2 q ,0,m 2 q ) - B l (m 2 q ,0,m 2 q )) 

+ sm29gmg(B (m 2 q ,ml,m 2 Sl ) - B (^,m|,m|)) 
+m q (Bi(m 2 q ,ml,m 2 q - 1 ) +5 1 (m 2 ,m?,m|))] , (21) 

which is derived from the graphs of Fig. Id. Furthermore, in the phase space term 
ftpoie in (20) we have to take the pole squark masses given by 



™!(p° le ) = m l(Q)m; ~ U U™>1) 

= <(Q)m + \^ [-3^(0) + 4m| B (ml , 0, m|) + 2m|5 1 (m| , 0, m 2 ) 



cos 2 26 q A{mi) - sin 2 20 f A(mf J + 4A( 



+4m 2 gB (m 2 qi ,m 2 g , m 2 ) + 4m|.Si(m|, mf, m 2 ) 

-(-) i - 1 4sin2^m g m §J B (m 2 i ,m 2 ,m 2 )] . (22) 

4 Gluon emission 



The infrared divergences in ( P0"D are cancelled by including the 0(a s ) contribution 
from real gluon emission from t and b (Fig. le). The decay width of H + (p) — > 



ti(ki) + bj(k 2 ) + g{kz) is given in terms of the DR parameters as 
T(H+ - Ubjg) = ^^^[( m ^- m |- m |)J 12 - TO ? J n - TO ? Jaa-Jx-Ja]. (23) 



The functions /„, and 7„ m are defined as [12 



1 f d 3 k 1 d 3 k 2 d 3 k 3 4 



~ W 2E X 2E 2 2E, 5 (P - *i - *2 - ^) (2 ^ i + A2) ___ (2 ^ + A2) - (24) 



The explicit forms of i^...^ are given in [I5|. In (J?3|), 111,22.12 are infrared divergent. 



We have checked that the infrared divergences in (|23|) cancel those in ((2^). In the nu- 
merical analysis we define the corrected decay width as Y corr (H + — > Ubj) = T(H + — > 

t^O + r^+^tJ^). 

5 Numerical results and conclusions 

As in ref. [0, we choose { m#, m fi &(pole), M, /x, tan/3, Mq, A } as the basic in- 
put parameters of the MSSM, taking M = (a 2 /a s )rrig = (3/5 tan 2 9 W )M', Mq = 
Mq(Q)vr = M {Q) m = M^(Q) m = M~ L and A = A t (Q) m = A b (Q) m = A T . Here 
M (M f ) is the SU(2) (U(l)) gaugino mass, a 2 = g 2 /4ii, and (Mg, A T ) are the mass 
matrix parameters of the slepton sector |J. The parameters M, M', Mg, and A T 
do not receive 0(a s ) QCD corrections. The theoretical and experimental constraints 
for the basic input parameters are described in ref.||. We take mz = 91.2GeV, 



m w = 80GeV, m t (pole) = 180GeV @, m 6 (pole) = 5GeV, sin 2 6 W = 0.23 and 



a 2 = a 2 {m z ) = a/sm 2 6 w = (l/129)/0.23 = 0.0337. For the running QCD coupling 
at the renormalization scale Q, a s = a s (Q), we always take the one-loop expression 
a s (Q) = 127r/{(33-2n / )ln(Q 2 /A 2 / )}, witha s (m z ) =0.12, and the number of quark 



flavors rif = 5(6) for nib < Q < m t (for Q > m t ). 

We define the QCD corrections as the difference between the 0(a s ) corrected 
width r^ orr = T COTT (H+ -> t-bj) of eqs.(20) plus (23) and the tree-level width T\f e = 
r(°)(#+ -> Ubj) of eq.(6) where (m,(pole), Mq, A) are substituted for (m q , Mgfjf), 
A q ). The QCD corrections depend on the renormalization scale Q. We choose the 
optimum value of Q (Q op t) such that A(Q) = X! ( m <?;(t ree ) — m gi(P°l e )) 2 is 

q=t,bi=l,2 

minimized, where m^(tree) refers to the mass defined by (1-4) calculated with 
m g (pole), Mq, and A. 

In order not to vary too many parameters, in the following we fix fi = 300GeV, 
and take the values of M and tan f3 such that m^o ~ 50GeV as in || where Xi 
is the lightest neutralino. In Fig. 2 we show the m# dependence of the tree- level 
and corrected widths T tree {tb) = ]T Tg ee and r corr (t1) = J2 T T> and the tree " 

i,i=l,2 _ i,i=l,2 

level branching ratio B tTec (tb) = J2 B tTCC (H + -> t-b 3 ) § for (a) Mq = 250GeV, 
A = 650GeV, tan/3 = 2, M = 120GeV, and (b) Mq = 136GeV, A = 260GeV, 
tan/3 — 12, M — HOGeV. In these two cases we have (in GeV units): (a) Q opt = 
216.5, m g = 380, {m h , m i% , m~ bi , m h ){ (tree), (Q opt ) m , (pole)} = {(60,429,251,254), 
(70,420,251,254), (58,428,258,262)}, m-+ = 94, and (b) Q opt = 213.5, m g = 349, 
K- 1 ,m f2) m Sl ,m g2 ){(tree) ) (Q opt ) 55 ,(pole)} = {(81,302,62,193), (74,292,120,164), 
(66, 302, 126, 171)}, m~+ = 98. Here xt is t ne lighter chargino. In both cases we see 
that the tb mode dominates the H + decay in a wide nin+ range at the tree level, and 
that the QCD corrections to the tb mode are significant, but that as a whole they do 
not invalidate the tb mode dominance. 
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In Table 1 we show the values of the B trcc (tb), the QCD corrections C = 
(T corr (tb) - T trcc (tb))/T trcc (tb), and Qj = (r^ orr - T\f c )/T\f c for typical values of Mq 
and A, for (a) m H+ = 400GeV, tan/3 = 2, M = 120GeV, and (b) = 400GeV, 
tan/3 — 12, M — HOGeV. We see again that the tb mode dominates in a wide region 
also when the QCD corrections are included. The QCD corrections can be very large 
at some points of (Mq,A); e.g. C = -0.692 at (175GeV, OGeV) and C 12 = 0.734 
at (225GeV, — 350GeV) in Table lb. This occurs when m# ~ uif. + m~ b ^. This 
enhancement is just a kinematical effect due to the QCD corrections to m 4 ~. and m^ . 

In conclusion, we have calculated the 0(a s ) QCD corrections to the decay 
width of H + — > tibj, including all quark mass terms and — qn mixing. We find that 
the QCD corrections are significant but that they do not invalidate our previous con- 
elusion at tree-level about the dominance of the tb mode in a wide MSSM parameter 
region. 
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Table 1 

B tme (tb), C, and C„ for typical values of Mq and A, for (m^(GeV), m t (pole)(GeV), 
M(GeV), /i(GeV), tan/3)= (400, 180, 120, 300, 2) (a) and (400, 180, 110, 300, 12) (b). 
The LEP bounds m- z -p^45GeV and the requirement m^^i > m^— 50GeV) are 
satisfied. Q opt denotes the optimum value of the renormalization scale Q as defined 
in the text. The * in the column of CV, means that the tjbj channel is not open at 
the tree level and/or one-loop level. In Table lb the columns of C 2 i and C22 are 
omitted, because the corresponding channels are not open, except C21 = —0.141 at 
(Mq,A) = (150GeV,0GeV). 
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(a) 





A(GeV) 


B tTee (tb) 


C 


Cii 

w 11 


(_y 1 o 

w 1Z 


w Z 1 


WOO 

w zz 


Q nnt (GeV) 

^ opt V / 


75 


o 


0.715 


0.054 


-0.251 


0.143 


0.025 


0.664 


224.0 


75 


200 


0.793 


0.158 


-0.090 


0.363 


-0.211 


0.319 


225.0 


75 


300 


0.818 


0.161 


-0.181 


0.261 


-0.229 


0.363 


240.5 


100 





0.690 


0.090 


-0.229 


0.165 


0.115 


1.051 


222.0 


100 


300 


0.788 


0.219 


-0.166 


0.295 


-0.181 


0.509 


234.5 


125 


100 


0.689 


0.240 


-0.170 


0.332 


0.127 


0.891 


217.5 


125 


350 


0.702 


0.175 


-0.220 


0.267 


* 


* 


238.0 


150 


300 


0.607 


0.280 


-0.124 


0.402 


* 


* 


221.0 


150 


400 


0.719 


0.159 


-0.272 


0.232 


* 


* 


238.0 


200 


500 


0.730 


0.188 


-0.375 


0.215 


* 


* 


226.5 


250 


650 


0.752 


0.209 


-0.437 


0.236 


* 


* 


216.5 



(b) 



M^GeV) 


A(GeV) 


B trec (tb) 


C 


C n 


C\2 


Qopt(GeV) 


150 


-200 


0.736 


-0.152 


-0.200 


0.531 


157.5 


150 





0.715 


0.062 


-0.383 


* 


159.0 


150 


250 


0.740 


-0.013 


0.047 


0.394 


200.0 


175 


-250 


0.651 


0.120 


-0.132 


0.399 


167.5 


175 





0.612 


-0.692 


-0.432 


* 


169.0 


175 


300 


0.681 


0.251 


0.139 


0.323 


204.5 


200 


-300 


0.606 


0.169 


-0.061 


0.395 


166.5 


200 


350 


0.656 


0.306 


0.239 


0.345 


201.0 


225 


-350 


0.512 


0.352 


0.012 


0.734 


161.5 


225 


400 


0.588 


0.501 


0.355 


0.587 


193.5 


250 


500 


0.616 


0.383 


0.579 


0.319 


193.5 
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Figure Captions 

Fig.l All diagrams relevant to the calculation of the 0(a s ) SUSY QCD corrections 
to the width of H + -> tjbj in the MSSM. 

Fig.2 The m H dependence of r tree (£ t^) = r tree (t6) (dashed line), r corr (£t~4) = 

r corr (t6) (solid line), and B tiee (tb) (short-dashed line) for (m t (pole)(GeV), M(GeV), 
Ai(GeV), tan/3, M^GeV), A(GeV))=(180, 120, 300, 2, 250, 650)(a) and (180, 
110, 300, 12, 136, 260)(b). T^iUbj) = T^ cc (dashed lines) and V^iUbj) = 
r? orr (solid lines) are separately shown for =(1,1) and (1,2). 
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Fig. 2 
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